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Abstract: While automatic detection of point sources in astronomical images has experienced a great 
degree of success, less effort has been directed towards the detection of extended and low-surface 
brightness features. At present, existing telescopes still rely on human expertise to reduce the raw data 
to usable images and then to analyse the images for non-point like objects. However, the next generation 
of radio telescopes will generate unprecedented volumes of data making manual data reduction and 
object extraction infeasible. Without developing new methods of automatic detection for extended 
and diffuse objects such as supernova remnants, bent-tailed galaxies, radio relics and halos, a wealth 
of scientifically important results will not be uncovered. In this paper we explore the response of the 
Circle Hough Transform to a representative sample of different extended circular or arc-like astronomical 
objects. We also examine the response of the Circle Hough Transform to input images containing noise 
alone and inputs including point sources. 
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1 Introduction 

The radio sky is rich in scientifically important ob- 
jects containing circular or arc-like structures. These 
objects range in physical size from the shells of dying 
stars to shocks created in merging clusters of galaxies. 
On the small scale, objects such as supernova rem- 
nants stereotypically contain complete or partial circu- 
lar she lls, or in cases when their sym metry is distorted, 
both rfMineshig e and Shibata|[T990t ). Arc- like features 
are also found in Active Galactic Nuclei ( AGN) , either 
in the overall s ource shape as in the case of bent-tail 
radio galaxies ([Johnston- Hollitt et a l. 2004; Mao et all 
2009), or as features observed within the lobes of the 
AGN themselves (jFeain et al.1 1201 lh . Finally, giant 
Mpc-scaled radio relics, which occur on the edges of 
merging clusters, are found as single or double arc-lik e 
regions (|Johnston-Hollittll2003l : iBonafede et al.ll2009l ). 
In addition to their arc-like morphology, these objects 
are often low in surface brightness and exhibit breaks 
or gaps which make them hard to detect in an auto- 
mated fashion. This is particularly true of techniques 
that rely on either high signal-to-noise or continuous 
edges in order to define the full extent of the source. 

In order to reliably detect broken, low-surface bright- 
ness sources we require algorithms which do not rely 
solely on the behaviour of individual pixels, but exploit 
additional non-local information. The advantages of 
such algorithms are that they are robust in the pres- 
ence of non-connectedness and noise. Examples of such 
techniques include template matching, wavelet and re- 
lated integral transforms and Hough transforms. As 



the Hough transform is known to be robust to the pres- 
ence of partia l or occluded target objects and to noise 
(|Davieslll988l ). it is a natural choice to explore for the 
automatic detection of circular features in astronomi- 
cal data. 

There has been some previous use of the Hough 
transform in astronomy. It has be en used to remove 
contamination from optical surveys (Storkev et al. 2004), 
to characterise echelle orders ([Ballesterl 1 19941, P.996I) , 
and to characterise craters on the moon (|Jahnlll994r ). 
However, it has not been previously applied for either 
source detection or characterisation in radio astron- 
omy. The increased resolution and field of view that 
will become routine when using instruments like the 
Square Kilometre Array (SKA) (|Schilizzi et al.ll2008T ) 
and its precursors such as the Australian SKA Pathfinder 
(ASKAP) and the Murchison Widefield Array (MWA), 
force a paradigm shift towards greater reliance on au- 
tomated data extraction techniques such as that de- 
scribed here. 

In this work we concentrate on the detection prob- 
lem; that of locating circular regions of radio emission 
that are unlikely to have arisen from statistical fluctu- 
ations in image noise. We also believe that the tech- 
niques described in this paper could contribute to the 
related characterization problem, where morphologi- 
cal features of the uncovered emission are accurately 
determined. As a result in this paper we have concen- 
trated on detecting circular structures, rather than the 
accuracy or precision with which they are described. 

In this paper we present a description of the cir- 
cle Hough transform (CHT) in Section [2] followed by 
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a detailed discussion of the applications of the CHT 
as a detection algorithm for use in Radio Astronomy 
in Section [3] Section [4] presents a discussion of com- 
putational issues and response to image artifacts and 
presents a suggested framework for the practical imple- 
mentation of the CHT for next generation radio tele- 
scopes. We present a summary of the work and final 
conclusions in Section [5] We have adopted a stan- 
dard set of cosmological parameters throughout with 
Ho 73 kms -1 Mpc"\ Q m 0.27 and Q A 0.73. 

2 Circle Hough Transforms 

The family of algorithms known collectively as Hough 
transforms are capable of finding and characterising a 
variet y of geometrical objects in image or image like- 
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data (jDuda and Hartlll972l ; lKerbvson and Athertonll995h . 
The application of the Hough transform for circle de- 
tection was a very early developm ent in the field, and 
has since found wide application (Kimme et al1 ll975l ; 
llllingworth and Kittled! 19881 : iRadet al.ll2003T ). 

In its simplest form the circle Hough transform is 
computationally challenging, with both computational 
effort and memory consumption scaling as 0(n 3 ) when 
transforming an n x n pixel image. Many variants of 
the CHT algorithms have consequently been developed 
to impro ve the compu t ational performance of the al- 
gorithm db avies 1988; Kerbvson and Athertonl Il995l : 
Ioannou et al. 1999; Chiu et al. 2010), but current im- 
plementations of the CHT remain confined mostly to 
offline processing. 

As is the case for all Hough transforms, the CHT 
seeks to map an image into a quantised parameter 
space that describes the target objects. An arbitrary 
circle in two dimensions can be described by 



(x - x ) 2 + (y- yof - r 2 



0, 
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where xo and yo are the x and y coordinates of the cen- 
tre of the circle and r is the circle's radius. Each differ- 
ent three-tuple x := (xo,yo,r) uniquely parametrises 
different circles that may be present in an image. Each 
point in the Hough space therefore corresponds to a 
complete circle with a particular location and size. The 
operation of the CHT is to map the input image into 
the parameter space such that the transform has a high 
value for parameter tuples consistent with objects that 
are present in the image. 

In its common form, the Hough transform oper- 
ates on each image pixel independently. For each pixel 
the algorithm finds the set of parameter tuples that 
would result in an object passing through the consid- 
ered pixel. That is, consideration of an image pixel 
f(xi,yi) finds the set V of all (xo,yo,r) tuples that 
satisfy equation [T] 

V = {(x ,2/o,r) : (xi - x ) 2 + (yi - y ) 2 - r 2 = 0} 

(2) 

Each location in Hough space that corresponds to 
a member of V is then incremented. Thus each image 
pixel can be regarded as casting a set of "votes", for 
those parameter tuples that are consistent with itself. 
Examination of equation [2] reveals that the set of votes 
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Figure 1: Formation of the Hough transform via 
convolution. The input image f(x,y) is extended 
to three dimensions and then convolved with the 
vote distribution K(x) to generate the Hough 
transform H(x). For clarity we have represented 
the Hough cone as smooth, though it is discretized 
on the same scale as the input image. 



arising from a single pixel forms a right cone i n Hough 
space (|Duda and Hartll 19721 : lYuen et al.lll990T ). 

The vote distributions from multiple image pixels 
are added and produces a high number of votes for pa- 
rameter combinations that describe geometric figures 
that are actually present in the image. The operation 
of the CHT can be conveniently regarded as the three- 
dimensional convolution of the input image with the 
conical vote distribution (Hollitt 2009), as illustrated 
in Figure [U 

For our application, if we calculate the CHT of 
a radio image and then find the peak in the Hough 
transform we will recover the radius and center of any 
circular astronomical features in the input radio map. 

2.1 Detection Procedure 

Peaks in the Hough transform of an arbitrary input 
image can arise even in the absence of a circular region 
of diffuse emission. In this section we briefly discuss 
the production of peaks due to noise in the input image 
and note that peaks are also expected because of the 
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presence of point sources. 

Let us first consider the expected Hough peak size 
that would arise in Hough space as a result of input 
image amplitude noise, which we expect to be nor- 
mally distributed with mean d and standard devia- 
tion s (that is the input image is ~ Af(a,s 2 )). When 
searching for a circle having radius r the Hough trans- 
form adds the magnitudes of approximately 2?rr pix- 
els. For a random input image we would therefore ex- 
pect that the distribution of magnitudes for the Hough 
transform would be ~ Af (27rrd, 27rrs 2 ) . 

One can therefore select an appropri ate t hresh- 
old for the source detector as 2nra + £sv / 27rr, for £ 
the number of standard deviations above the expected 
mean that one wishes to use. For an image a few 
hundred pixels in extent, f = 5 should ensure that 
Hough peaks arising solely from noise are rarely en- 
countered. A considerably more aggressive threshold 
could be used if a higher false- positive detection rate 
were acceptable for a particular application. 

The noise in radio images arising from surveys is 
unlikely to be stationary, so the appropria te values of 
a and s will need to be calculated locally (|Huvnh et al.l 
2011). Note also that the analysis above would need 
to change slightly if the Hough transform were used 
to search for structures in polarimetric surveys, as the 
noise in polarisation images is Rician rather than Gaus- 
sian. 

We would also expect that the Hough transform of 
an arbitrary radio image will contain a number of pu- 
tative circles due to the presence of point sources. In a 
noiseless image the Hough transform would find circles 
centered half way between each pair of point sources. 
Such circles will have a diameter approximately equal 
to the separation of the points sources. However in a 
noisy image the circle is typically slightly displaced to 
include high noise pixels in its circumference, as well 
as the two point sources. 

A circle can also be fitted though any set of three 
points (with the exception of three collinear points), 
so we would also expect that the Hough transform of 
an arbitrary field would contain a peak corresponding 
to each combinations of three point sources. However, 
in practice many such peaks would correspond to very 
large circles, or circles with centres lying well outside 
the image boundaries. That is, if the image contains 
N point sources then the Hough space will include 
peaks. In the case of perfect point sources, the magni- 
tude of these peaks would be the sum of the amplitudes 
of the three selected sources. In practice they will de- 
viate somewhat from this due to the non-ideal beam 
pattern. 

As the magnitudes of the Hough peaks arising from 
point sources are dependent on the source intensities, 
it is impossible to remove their effects with a simple 
threshold. Instead it will be necessary to examine cir- 
cles detected with the Hough transform to determine 
whether they have arisen via this mechanism. We have 
not yet examined this problem in detail, though it will 
be an important component of any system intended 
to automatically process survey data using the Hough 
transform. 

Combinations of more point sources are of course 



also possible if they happen to lie on a common cir- 
cle, though this becomes increasingly unlikely as the 
number of points increases. 

The overall procedure proposed to find peaks in 
Hough space is as follows: 

1. Find the peak Hough transform amplitude for 
each value of r. 

2. Calculate how many standard deviations above 
the expected noise-only level 2irra the maximum 
value of the Hough transform lies for each r. 

3. Search for local maxima in the resulting ampli- 
tude as r varies. Such local maxima correspond 
to circles in the input image. 

4. Check for circles that obtain most of their votes 
from two (or more) point sources and remove 
them from the list of possible diffuse emission 
candidates. 

An example of this procedure is provided in Section 

3 Results 

We applied the CHT to images containing a variety of 
arc-like radio sources in order to test features of the al- 
gorithm including filtering and robustness in the pres- 
en ce of no i se. Th e variant of the CHT used is described 
in iHollittl (|2009h , though the results described below 
are not dependent on the particular implementation 
of the algorithm. Each of the three species of source 
tested — SNRs, tailed radio galaxies and giant radio 
relics — represent populations which will be present in 
next generation radio surveys. It is expected that im- 
portant science will be determined from the statistical 
characterisation of such populations, making reliable 
automatic classification an important goal. 

3.1 Response to non-diffuse struc- 
tures 

We expect that the distribution of the Hough trans- 
form of a Gaussian noise field - A/"(a,s 2 ) will be cen- 
tered on 2nra with standard deviation s\f2nr. This 
allows us to determine an appropriate threshold to as- 
sess the reliability of any proposed circle detections. 

Gaussian noise fields of 243 x 243 pixels were pre- 
pared with the pixels values distributed with known 
mean and variance. The Hough transform was then 
calculated to find the distribution of the magnitudes 
in Hough space. For each value of r the largest value 
in the Hough transform was found and compared with 
the theoretical distribution. 

Fig [2] shows the results of one such simulation, hav- 
ing the input pixels distributed as ~ J\f (1,0.25). As 
can be seen the largest values for the Hough trans- 
forms typically lie between three and four standard 
deviations from the expected value of the transform. 
This is consistent with the expected extreme value of 
243 2 samples drawn from a normal distribution. 

We therefore conclude that for the image size of 
243 2 pixels used throughout this paper, a threshold of 
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Figure 2: The red crosses indicate the largest value 
of the Hough transform, taken across the 243 2 val- 
ues for each target circle radius. The bold blue line 
indicates the expected value for the Hough trans- 
form of the random input field, with the thinner 
lines indicating progressive steps of one standard 
deviation in the expected value of the Hough trans- 
form. 



five times the input image standard deviation is suffi- 
cient to exclude features in Hough space arising from 
image noise. 

Figure [3] shows the result of applying the Circle 
Hough Transform to an im age from the Molon glo Galac- 
tic Plane Survey (MGPS) (|Green et alJI 19991 ) contain- 
ing no known diffuse structures, but a number of iso- 
lated point sources. As discussed in 12.11 the Hough 
transform results in a peak corresponding to a cir- 
cle centered roughly between the two brightest point 
sources. It also finds a second circle that additionally 
passes through a third point source. 




Figure 3: Response of the Circle Hough Trans- 
form to a collection of point sources. The centre 
of the solid circle lies between the two bright point 
sources. The dashed circle also passes through a 
third point source to the left of the image. 



3.2 Supernova Remnants 

When a massive star at the end of its life becomes a 
supernova, an energetic shock wave is released from 
its surface. As the shock travels outward it acceler- 
ates electrons in the interstellar medium to produce 
emission at radio frequencies. In the simplest case 
the resulting emission would be spherically symmet- 
ric as the spherical shock passes through a homoge- 
neous interstellar medium. However, where the shock 
wave is asymmetric, or when there is exiting struc- 
ture in the gas cloud surrounding the supernova, then 
more complicated supernova remnants are pr oduced 
(|Gaensler and Slandl2003 : lOrlando et "aDl2007h . 

Our galaxy contains a large number of known su- 
pernova remnants (SNRs); the result of supernovae oc- 
curring in roughly the last million years. The Molonglo 
Observatory Synthesi s Telescope Supern ova Remnant 
Catalogue (MSC) (|Whiteoak and Gr een 1996) contains 
a list of 75 SNRs from the 843 MHz survey of the 
Southern Galactic plane within the area 245° < / < 
355° in Galactic longitude and b = ±1.5° from the 
Galactic plane. Statistical studies of the SNR pop- 
ulation allows much more information to be inferred 
than would be possible by examination of one SNR 
alone. In addition to the valuable information gleaned 
on the average nature of the ISM surrounding super- 
novae, statistical analysis shows that there is presently 
a significant discrepancy between the total number of 
know n SNRs and that expected from stellar lifecy- 
cles ([Brogan et al.l 12006). Comparison with the ex- 
pected rate of supernovae production from OB star 
counts, pulsar creation rates, iron abundance and ob- 
servation of supernova in other galaxies suggests that 
there sho uld be over 1000 SNRs detectable in th e Galac- 
tic plane ([Li et al.lll99ll : iTammann et al.lll994T ). How- 
ever, at present only ^250 have been found despite 
an increasing number of surveys of the region (|Greenl 
120041 : IBrogan et aD 120061 ) . The detection of the pop- 
ulation of SNRs is therefore an important aspect of 
radio surveys of the Galactic disk. 

The missing SNRs are believed to be those that are 
young and distant and those whic h are old and nearby 
([Greenl |2004| ; IBrogan et al.l 12006) . In the former case 
the issue is simply a lack of resolution in current sur- 
veys of the Galactic plane leading to such small an- 
gular sized, bright sources being misclassified as point 
sources. This issue will be readily addressed and re- 
solved by the next generation of radio telescopes such 
as the SKA. In the latter case however, the SNRs will 

be both faint and large and if the prese nt manual "eye- 

balling" techniques continue to be use d (White oak and Greenl 
119961 : iFilipovic et~aD 120021 : IBrogan et al.l 120061 ) thev 
are likely to remain elusive. Additionally, manual in- 
spection is error prone and often re sults in objects be- 
ing missed in the in itial data release (Whit eoak and Greenl 
Il996l ; iGreenl l2004f ) . Finally, while this approach has 
been adequate (although not ideal) in the past, fu- 
ture generations of radio telescopes will generate such 
a vast stream of data that manual inspection will not 
be practical. 

We have selected a subset of four representative su- 
pernova remnants from the MCS catalogue, as shown 
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Figure 5: Slices through the Hough transform 
of G337.3+1.0. Slices are taken through several 
planes corresponding to a set of different possible 
SNR radii, r = (16,32,48,64,80) pixels. 



in Fig 21 This subset was selected to range from a 
near complete ring J4ji) , through two partial rings JUd 
and [4]:) to a distorted example J4jl). These four SNRs 
were intended to provide an increasing challenge for 
the algorithm. 

A CHT algorithm was performed on each of the 
fou r supernova r emnants using the procedure outlined 
in ([Ho llitt 2009). The four input data sets were each 
padded to 243 pixels square. A Hough transform cal- 
culation is typically preceded by an edge detection step 
to highlight the boundaries of objects. However, su- 
pernova remnants are inherently shell-like, so we found 
the edge-detection step to be unnecessary. Compari- 
son of the transformation of the raw and edge-detected 
data resulted in differences in only a few pixels in the 
estimated SNR position and radius. 

Figure [5] shows the Hough transform of one of the 
SNRs. For clarity we selected the almost complete ring 
of G337.3+1.0. While the Hough transform results in 
a full three-dimensional array this is difficult to present 
graphically. The figure therefore shows the magnitude 
of the Hough transform obtained for a small subset 
of the possible SNR radii. As can be seen for small 
radii the intensity pattern is well distributed, indicat- 
ing that there is no strong peak for these radii and 
hence no small radius circles in the image. However, 
in the r = 64 plane we can see a strong peak toward 
the centre of the field. This indicates that the input 
image contains a circle with radius approximately 64, 
which is located toward the centre of the image. For 
larger radii we again see the intensity pattern is spread, 
so there are no larger circles present. 

The Hough transforms of the four SNRs were searched 
to find their maxima. As an example of the procedure 
used, Fig [6] shows the peak value of the Hough trans- 
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Figure 6: Maximum values of the Hough trans- 
form of G337.3+1.0 as a function of radius are indi- 
cated by the red crosses. The thick lower line indi- 
cates the Hough transform expected from the im- 
age noise alone, with the increasingly higher lines 
indicating one standard deviation in the expected 
noise-only values. 

form of G337.3+1.0 as a function of the target circle 
radius. Also shown in the figure are the peak values ex- 
pected if the image contained noise alone. The figure 
clearly shows that the input image contains features 
other than noise. 

The height of the peak above the noise-only level 
is then calculated in units of standard deviations, that 
is, we calculate the statistical significance of the Hough 
transform. The result of this calculation is shown in 
Figure Local maxima in this graph indicate peaks 
in Hough space, which can then be used to extract the 
location of circles in the original image. In the case 
of Figure the prominent peak at a radius of 65 pix- 
els corresponds to the shell of the supernova remnant. 
The local maxima around 10 pixels is associated with 
the circular feature in the lower right portion of the 
SNR, as seen in Figure |4Jl. 

These peaks in the transforms correspond to the 
circle that is most consistent with the input data in 
each case. Extraction of the coordinates of the peaks 
therefore results in an estimate of the location and 
radius of the supernova remnants. The effectiveness of 
the CHT in identifying the circular structures in the 
four supernova remnants can be seen in Figure HJ In 
all four cases the best fitting circles returned by the 
algorithm are in good agreement with the radio data. 

3.3 Bent-tail Radio Galaxies 

Bent-Tail (BT) galaxies are Fanar off- Riley class I /II 
(jFanaroff and Rilevlll974 ) radio sources that consist of 
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Figure 4: Four supernova remnants taken from the MCS catalogue ([Whiteoak and Greenl fl996). On the 
left panel subfigure a. is G337.3+1.0, b. is G302.3+0.7, c. is G315.4-2.3 and d. is G317.3-0.2. On 
the right panel are the same four objects overlaid with circles that indicate the features detected via the 
Hough Transform of each object. 
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Figure 7: Maximum value of the Hough transform 
of G337.3+1.0 as a function of radius. The ampli- 
tude is expressed as the number of standard devi- 
ations between the peak value and the value ex- 
pected from a field consisting solely of Gaussian 
noise. 



a core, coincident with the optical host galaxy, and two 
tails of bent radio emission. The bending of the radio 
jets i s believed to resu lt from the actions of ram pres- 
sure (|Gunn and Gottlll972t ) due to the relative motion 
of the galaxy and its surrounding. In the classic case 
this results in a 'C ' shaped radio object such as SUMSS 
J032752-532613 ([Johnston-Hollitt et al.ll2004T ) seen on 
the left of Figure [8j though more complicated struc- 
tures such as PKS J0327— 532 seen on the right of the 
same figure do occur. 

BTs have been found in the densest and most tur- 
bulent regi ons of galaxy clu sters in the local (z < 0.07) 
Universe ([Mao et al.ll2009f ) at a rate of about 1-2 per 
cluster. BTs are thought to act as signposts for over- 
densities in large-scale structure, with associations be- 
tween tailed- ga laxies and clusters pers isting out to at 
least z rsj 2 (Dehgh an et al.ll201lL I2012T ) making them 
an important class of objects for science with future 
all- sky radio survey. 

One of the most important parameters to extract 
when using BT galaxies to obtain physical informa- 
tion on the environment in which they are embedded 
is the radius of curva ture of the source ([Freeland et al.l 
120081 : lFreelan3l2010T ). This technique is only applied 
to wide-angle tailed sources (WAT) such as SUMSS 
J032752— 532613, so it is typically performed subse- 
quently to the initial source detection phase. For the 
test described in this section we wished only to find the 
radius of the WAT, rather than extract information 
about the many wiggles in the tails of the other BT. 
Consequently we imposed a constraint on the allowable 
circle radii before running the transform. This demon- 
strates the ability of the CHT to search for sources of 
only a certain angular size and to perform the simul- 
taneous science-driven detection and charactisation of 
the WATs. The successful results are seen in the left 
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Figure 8: 1.4 GHz image of a pair of head- 
tailed radio galaxies located in the centre of Abell 
3125 (iJohnston-Hollitt et all 120041 ). The classic 
'C shaped HT (SUMSS J032752-532613) has been 
correctly identified by the CHT over all of the 
other radio sources and artifacts in the image. 

of Figure [8] 

As discussed in more detail in Section 14.21 radio 
telescopes with only a few elements suffer from gaps 
in the uv-plane which result in circular artifacts in the 
image-plane. In the case of the im a ge us ed here, the 
image, from I Johnston-Hollitt et~al . ( 2004), comprises 
of a mosaic of 4 pointings of 1.4 GHz data from 1.5 and 
6 km configurations of the pre- C ABB Australia Tele- 
scope Compact Array (ATCA) which is known to have 
significant gaps in the distribution of baseline lengths. 
The size of the resultant artifacts in the image plane is 
a function of the antenna element spacing and is there- 
fore known in advance allowing us to employ filtering 
to prevent the CHT returning artifacts as sources. As 
next generation radio telescopes will employ a larger 
number of elements, this is likely not to be a major 
issue for future surveys. 

3.4 Radio Relics 

Radio relics are a rare class of large arc- like radio 
sources seen on the periphery of galaxy clusters. They 
are believed to be the result of Shockwaves which oc- 
cur when clusters of galaxies merge, the most powerful 
events since the Big Bang. The physical size of these 
sources can be up to six million light year across, such 
as in the northern relic in Abell 3667 seen in Figure 
[9] Because of their extremely low surface brightness 
they are very difficult to detect with conventional pixel 
threshold based searches. As in the case of faint SNRs, 
radio relics should be more easily detected with al- 
gorithms utilising non-local information such as the 
CHT. 

We tested the CHT on th e field of Abell 3667 us - 
ing 1.4 GHz ATCA data from lJohnston-Hollittl (|2003T ). 
These data use five different configurations of the tele- 
scope and so have good sampling of the uv-plane and 
hence a minimum of circular artifacts. This particular 
image has had a taper applied in the uv-plane to give a 
resolution of approximately one arcminute. The data 
contain a number of extended sources including the 
double relics for A3667 that are located to the top right 
and bottom left of the image, two BT galaxies in the 
centre of the field and at the lower left above the relic 
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Figure 9: 1.4 GHz Australia Telescope Compact 
Array (ATCA) mos aic image of the doubl e relic 
cluster Abell 3667 ([Johnston- Hollittl l2QQ3h . this 
image covers approximately four square degrees on 
the sky and shows an example of one of the largest 
arc-like sources known. The CHT has identified 
the arc- like shape of the northern relic as well as bi- 
furcations in the intensity structure of the source. 



and a number of unresolved point sources. The CHT 
was employed without the filtering described in Sec- 
tion 13.31 and the resultant detected features in Hough 
space are displayed on the right side of Figure [9] 

Use of the CHT should isolate the location and 
scale of all circular features in a set of radio data to 
provide a catalogue of possible sources that can be as- 
sessed by an astronomer. In this case the algorithm 
correctly identifies the curvature of the northern relic 
and gives the approximate cluster centre (largest cir- 
cle). It also correctly characterises an internal curved 
structure in this relic (smallest circle). However, the 
characterisation of the relics is imperfect as Abell 3667 
has an elliptical, rather than perfectly circular struc- 
ture. Also the middle circle is influenced by the bright 
BT galaxy just north of the cluster centre which is 
unrelated to the relic. 



3.5 Robustness to Noise 

To illustrate the robustness of the CHT to noise we 
contaminated the supernova remnants in Figure [4] with 
increasing amounts of Gaussian noise and processed 
the resulting images until the SNRs could no longer 
reliably be detected. 

Figure [10] shows the supernova in Figure Ht con- 
taminated with additional noise until the Hough trans- 
form detector was 50% reliable. For this object the 
magnitude of the noise has been increased by a factor 
of 10 from that in the original image. As shown in 
the figure, the estimated location of the SNR in the 
contaminated image is slightly different to that in the 
uncontaminated image. However, as we are primar- 
ily interested in the detection of supernova remnants 
rather than their characterisation, we do not regard 
this as a significant problem. As might be expected 
the error in the suggested location of the object typ- 
ically decreases as the signal-to-noise ratio of the im- 
ages improves. 
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Figure 10: The input and results when the noise 
level in the supernova remnant from Figure Ht is 
increased until the Hough transform method is 
50% reliable. The solid line indicates the location 
of the SNR as estimated by the Hough transform. 
The dashed circle indicates the location of the SNR 
as predicted without additional noise added to the 
image before processing. 

4 Discussion 

4.1 Computational Issues 

Unfortunately the CHT is a computationally intensive 
algorithm, particularly when used to character ise cir- 
cles with unconstrained radii (|Yuen et alj fl990). The 
CHT is therefore likely to prove too slow to directly 
process the vast amounts of raw data that will be pro- 
duced by the next generation of radio telescopes such 
as the SKA. 

Instead, we propose a two tiered approach to the 
detection of circular and arc-like sources. The first 
stage will require the detection of regions of extended 
radio emission (that is emission from any objects other 
than point sources). This initial processing must be 
fast, but doesn't need to be very accurate. Some vari- 
ety of blob-detection algorithms would be suitable for 
this purpose. A fair amount of work has been carried 
out by the astronomical community o n this problem 
and s tandard algorith ms such as VS AD ([Condon et al.l 
119981 ) and Diichamp (|Whitindl20"l2h exist. Ideally the 
blob detection algorithm would additionally provide 
an estimate of the scale of a possible source. However, 
existing astronomical application of blob detection are 
poor at detecting broken large scale structures, so are 
likely to provide a poor estimate of both the size and 
centre of most SNRs or double relics. 

4.2 Artifacts 

Existing telescopes currently operating with v < 3 GHz 
necessarily use a relatively sparse sampling of spatial 
frequencies, which can result in circular calibration ar- 
tifacts in their output data. Such artifacts would con- 
found our detection algorithm in its simplest incarna- 
tion. However, as demonstrated in Section l331 grating 
lobes of known size can be filtered out before the fi- 
nal catalogue is returned. Additionally, many artifacts 
have circular features combined with radial spikes, al- 
lowing a combined CHT and linear Hough transform 
to be used to find and reject such sources. This would 
have the additional benefit of being able to identify 
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flawed regions in an automatically generated image be- 
fore a survey was released. 

5 Conclusions 

We have examined the Circle Hough Transform as a 
detection and characterization method for finding cir- 
cular or arc-like sources in astronomical data. We have 
tested the algorithm on MGPS and ATCA images of a 
range of sources including supernova remnants, tailed 
radio galaxies and radio relics which all represent scien- 
tifically important classes of radio sources for the next 
generation of radio telescopes. We have shown that 
the technique is successful at recovering such sources 
and providing physical information for them in an au- 
tomated way. 

At presen t thresholding algorithms such as Duchamp 
(|Whitind 120121 ) are able to produce catalogues of ra- 
dio sources based on the requirement of having a cer- 
tain number of contiguous pixels over a certain signal- 
to-noise level. These algorithms will usually return 
either the position of the brightest pixel or the geo- 
metrical centre of the source as the central coordinate. 
This approach will work adequately for the majority of 
symmetrical, unbroken sources. However, for arc- like 
sources such techniques will not produce the correct 
coordinates making automatic cross-matching to other 
wavelengths difficult (as is currently the case). Addi- 
tionally, with broken sources such as SNRs and double 
relics, such techniques will report multiple, unrelated 
diffuse sources (all with the wrong central coordinates) 
rather than identify a single source. 

Because it uses non-local information, the CHT is 
able to work with broken sources and is capable of 
detecting sources that have much lower signal-to-noise 
ratios than other algorithms. The CHT can therefore 
be expected to do a superior job of both detecting 
and characterizing some classes of sources, particularly 
those that are very diffuse. 

However, as the CHT is computationally expen- 
sive, it is unlikely to be sensible to run the algorithm 
over the entire sky for an all-sky survey such as EMU 
(|Norris et al.ll201lh . We propose that it would be bet- 
ter to use a conventional "blob detection algorithm" to 
preprocess the data before and then apply the CHT in 
certain parts of the sky only. These regions could ei- 
ther be selected based on criteria from the results of the 
conventional blob detection, or as regions where such 
low-surface brightness arc-like sources are expected such 
as in the Galactic plane or around galaxy clusters. 
This tiered approach is likely to be the most efficient 
method for generating a robust catalogue of diffuse 
radio sources and will be the subject of further exper- 
imentation. 

There are several improvements that could be made 
to the current CHT. Firstly, the robustness of the 
method can be improved by the development of a more 
effective method for searching the Hough transform for 
peaks. Simply finding the maximum value as we have 
done here is unlikely to be the optimal technique. Sec- 
ondly, as seen in some SNRs and double relics, the 
objects in question are elliptical rather than circular. 
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The Hough transform could be extended from the cur- 
rent circular implementation to an elliptical version. 
However this increases the number of parameters re- 
quired and hence would result in a five dimensional 
Hough space which would considerably increase the 
computation complexity. 

As has been noted, the Circle Hough Transform 
produces apparent circles as a result of point sources. 
If the Hough transform is to prove effective at auto- 
matic diffuse source detection, it will need to be com- 
bined with additional processing to remove the effects 
of point sources. Once this has been done the effective- 
ness of the combines detection system can be assessed 
by using it on a catalogue containing known diffuse 
structures. 

A catalogue of known supernova remnants has been 
manually compiled from th e second Molonglo Galac- 
tic Plane Survey (MGPS2) iMurphv et all (|2007T ). al- 
lowing ready comparison with the results of the auto- 
matic detection. Use of the CHT on such a dataset 
will allow the determination of the false negative and 
false- positive detection of the technique. This com- 
parison also admits the possibility of the detection of 
as-yet undiscovered supernova remnants and in par- 
ticular those in the class of being old and faint which 
have been predicted via several independent indicators 
but remain, as yet, undetected. Thus, in addition to 
characterising the false-positive and false- negative er- 
ror rates, performing the CHT on the MGPS2 could 
address one of the current problems in understanding 
the so-called "Galactic Ecology" . 
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